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Abstract 

In cyber-physical systems where sensors measure the temporal evolution of a given phenomenon of 
interest and radio communication takes place over short distances, the energy spent for source acquisition 
and compression may be comparable with that used for transmission. Additionally, in order to avoid 
limited lifetime issues, sensors may be powered via energy harvesting and thus collect all the energy 
they need from the environment. This work addresses the problem of energy allocation over source 
acquisition/compression and transmission for energy-harvesting sensors. At first, focusing on a single- 
sensor, energy management policies are identified that guarantee a maximal average distortion while at 
the same time ensuring the stability of the queue connecting source and channel encoders. It is shown 
that the identified class of policies is optimal in the sense that it stabilizes the queue whenever this is 
feasible by any other technique that satisfies the same average distortion constraint. Moreover, this class 
of policies performs an independent resource optimization for the source and channel encoders. Analog 
transmission techniques as well as suboptimal strategies that do not use the energy buffer (battery) or use 
it only for adapting either source or channel encoder energy allocation are also studied for performance 
comparison. The problem of optimizing the desired trade-off between average distortion and delay is 
then formulated and solved via dynamic programming tools. Finally, a system with multiple sensors is 
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considered and time-division scheduling strategies are derived that are able to maintain the stability of 
all data queues and to meet the average distortion constraints at all sensors whenever it is feasible. 

I. Introduction 

In the "smart world", wireless sensor networks (WSNs) play a central role in bridging the real 
and the digital worlds [DQ. WSNs are typically designed under the assumptions that communi- 
cation resources are limited by the energy available in the battery and that the most significant 
source of energy expenditure is radio transmission. However, modern cyber-physical systems are 
expected to operate over a virtually infinite lifetime. This can only be achieved by overcoming the 
limitations of battery-powered sensors and allowing the sensors to harvest the energy needed for 
their operation from the environment, e.g., in the form of solar, vibrational or radio energy [0, []3]. 
The regime of operation in which the system operates in a fully self-powered fashion is referred 
to as energy neutral [4]. Moreover, when sensors are tasked with acquiring complex measures, 
such as long time sequences of given phenomena of interest, and when transmission takes place 
over small distances, the energy cost of running the source acquisition system (sensing, sampling, 
compression) may be comparable with that of radio transmission flU, (fill- 
Based on the discussion above, in this paper, we address the problem of energy management 
for a WSN in which sensors are powered via energy harvesting and in which source acquisition 
and radio transmission have comparable energy requirements. We first focus on a system with a 
single sensor communicating to a single receiver, as shown in Fig. [TJ in order to concentrate on 
the main aspects of the problem. The sensor is equipped with a battery in which the harvested 
energy is stored. In each time slot, the sensor acquires a time sequence for the phenomenon of 
interest, which is characterized by a measurement signal-to-noise ratio (SNR) and autocorrelation, 
and stores the resulting bits, after possible compression, into a data queue. At the same time, 
it transmits a number of bits from the data queue to the fusion center over a fading channel 
with an instantaneous channel SNR. Based on the statistics of the energy harvesting process, 
and based on the current states of the measurement quality, of channel SNR, and of the data 
queue, the energy management unit must perform energy allocation between source acquisition 
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and data transmission so as to optimally balance competing requirements such as distortion of 
the reconstruction at the receiver, queue stability and delay. This optimization problem is the 
main subject of this work. We further extend our analysis to the problem of scheduling multiple 
sensors that are communicating to the same receiver. 

The model at hand is inspired by the work in [O, [|6) and 07]. In Q, the energy-harvesting 
sensor allocates power to data transmission over different channel SNRs, since the bit arrival 
process is assumed to be given and not subject to optimization. This is unlike our work in which 
a key problem is that of allocating resources between transmission and source compression in 
order to guarantee given constraints such as distortion and queue stability. The problem of energy 
allocation between source compression and transmission was instead first studied in [Q, fl6l, but 
in power- limited systems with no energy-harvesting capabilities. 

Other related works pertain to the study of energy-harvesting WSNs. This is a growing field 
with recent significant contributions. Here we only point to the works that are most related to 
ours, besides the ones already mentioned above. An information-theoretic analysis of a single- 
sensor system with energy-harvesting is presented in [§?], [HI, where it is shown that energy- 
harvesting does not affect the capacity of the channel, as long as one assumes that the battery 
has an arbitrarily large storage capacity. An optimal strategy for a single-sensor system that can 
control both the "acceptance rate" of the arriving bits and the power allocation with the aim of 
maximizing the throughput under stability constraints is developed in IfTOl . Optimal scheduling 
is instead studied in IfTTTl . |[T2Tl . |fT3l . The effect of a finite battery is studied in [[141 . where the 
trade-off between achievable rate and battery discharge probability is characterized. It is noted 
that all these works do not model the aspect of source acquisition and processing. 

The main contributions of this work are summarized as follows. (?) We propose a simple, 
but general, model for an energy-harvesting sensor operating over a time-varying channel (Sec. 
HU). (ii) For a single-sensor system, we design a novel class of distortion-optimal energy-neutral 
resource allocation policies that are able to stabilize the data queue and, simultaneously, to meet 
an average distortion constraint, whenever it is feasible by any policy (Sec. IIII-AI) . For the case 
where multiple sensors access the same uplink channel in time division, we identify a distortion- 



January 18, 2013 



DRAFT 



optimal energy-neutral class of scheduling policies (Sec. [V). (Hi) We compare the performance 
of the optimal policies with a number of less complex strategies, such as "analog" techniques 
lfl"5l (Sec. IIII-CI) and fixed time division multiple access (TDMA) scheduling strategies (Sec. 
IV-AI) . (iv) Finally, we formulate the problem of optimizing a desired trade-off between average 
delay and distortion, which is solved via dynamic programming tools (Sec. HVl) . 

II. System Model 

In this section, we introduce the system model, main assumptions and problem definition. 

We consider a system in which a single sensor communicates with a single receiver as depicted 
in Fig. [U The extension of this system to the case of multiple sensors and a single receiver is 
studied in Sec. [V] In most of the paper, we assume that the sensor performs separate source and 
channel coding, as described in the following. A different approach is considered in Sec. IIII-CI 

Time is slotted. The energy E k e R+ harvested in time-slot k is stored in an "energy buffer", 
also referred to as battery, with infinite size. For convenience, the energy E k is normalized to 
the number N of channel discrete-time symbols available for communication in each time slot, 
also referred to as channel uses. The energy arrival E k is assumed to be a stationary ergodic 
process. The probability density function (pdf) of E k is pe (e). The energy E k +i available for 
use at slot k + 1 is the residual energy from the previous slot plus the energy arrival at time-slot 
k + 1. This evolves as 



E, 



k+l 



E k — {T s ,k + T t 



t,k) 



+ E k+1 , (1) 



where T s k and T t k account for the energy spent in slot k per channel use for source acquisition 
and data transmission, respectively, as discussed below. Notice that the energy arriving at time 
slot k + 1 is immediately available for use in that slot. 

The sensor measures M samples of a given source during each slot. The quality of such 
observation in slot k depends on a parameter Q k e Q, which is assumed to be a stationary 
ergodic process over the time slots k. For instance, the sensor may perform measurements of 
the phenomenon of interest whose SNR Q k changes across blocks k due to source movement 
or environmental factors affecting the measurement quality. The set Q is assumed to be discrete 
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and finite, and the (stationary) probability mass function (pmf) for Q k is given by Pr(g) = 
Pr(Qfc = q), for q E Q. The sensor acquires the source in a lossy fashion. The loss, due to 
sampling, analog-to-digital conversion and compression, is characterized by distortion D k E M + , 
as measured with respect to some distortion metric such as the mean square error (MSE). 

The number of bits generated by the source encoder at the sensor at slot k is X k = f(Dk, T Stk , Qk), 
where / is a given function of the distortion level D k , of the energy per channel use allocated to 
the source encoder T s ^ and on the observation state Qk- The resulting bit stream is buffered in a 
first-input-first-output (FIFO) data queue with queue length X k . The function f(D k ,T Sjk ,Q k ) is 
assumed to be separately continuous convex and non-increasing in D k and T s k . For simplicity, 
we will denote such functions also as f q (D k ,T S:k ) = f(D k ,T Stk ,Q k = q). Some examples for 
function / will be provided below in Sec. III-AI 

The fading channel between sensor and destination is characterized by a process H k , assumed 
to be stationary ergodic, where H k E H, with set H being discrete and finite in order to ease the 
numerical evaluations. We assume a slowly time-variant scenario. The pmf of H k is given by 
Pr(/i) = Pr(H k = h), for h E K. The channel encoder uses the channel N times per slot, and 
the transmission requires T t)k energy per channel use. A maximum number g (H k ,T tjk ) of bits 
per slot can be delivered successfully to the destination. The channel rate function g (H k ,T tjk ) 
is assumed to be continuous, concave, non-decreasing in T t)k , and g(H k ,0) = 0. We also use 
the notation g h (T tjk ) = g (H k = h,T t k ). An example is the Shannon capacity on the complex 
additive white Gaussian noise (AWGN) channel g h (T tjk ) = N x log(l + hT t)k ) [16]. We remark 
that adopting the Shannon capacity implies the use of rate-adaptive schemes with sufficiently 
long codewords so that the block error probability becomes negligible. For the given function 
g(H k ,T t>k ), the channel encoder takes min[X k ,g (H k ,T t>k )] bits from the data buffer, using the 
selected transmission energy T t>k . Note that we do not consider the effects of channel errors nor 
the costs of channel encoding/decoding and of channel state information feedback, which are 
beyond the scope of the present paper and subject to future work. 
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Based on the discussion above, the data queue evolves as 



Xk — g (H k ,T tjk ) + f (Dk,T Sy k,Qk) ■ 



n + 



(2) 



To illustrate the trade-offs involved in the energy allocation between T tjk and T s>k , we remark 
that, by providing more energy T sk to the source encoder, one is able, for the same distortion 
level D k , to reduce the number / (D k , T s>k , Q k ) of bits to be stored the data buffer. At the same 
time, less energy T t)k is left for transmission, so that the data buffer is emptied at a lower rate 
g (H k ,T t k ). Viceversa, one could use less energy to the source encoder, thus producing more 
bits / (D k , T s>k , Qk), so that more energy would be available to empty the data buffer. 

A. Rate-Distortion-Energy Trade-Off 

In the following, we present some examples for function f q (D fc ,T S;fc ), as available in the 
literature. Recall that this function provides the trade-off between the distortion D k , the energy 
consumption T s fc and the number of bits produced by the source encoder. 



Example 1. Consider the observation model R kji = \fQ k U k ^ + Z kii , where M samples of the 

random process U k ,i, for i G {1, . . . , M}, are measured during the slot k and each measurement 
R k) i is affected by Additive White Gaussian Noise (AWGN) Z k with unitary variance. Parameter 
Q k represents the observation SNR in slot k. From [0, an approximated and analytically tractable 
model for /« (D k , T„, k ) is 



where b = N/M is the bandwidth ratio and f 2 (T 8tk ) = (x m&x[(bT Stk /T™ ax )~ 1/v , 1] models the 
rate-energy trade-off at the source encoder. The parameter £ > 1 is related to the efficiency of 
the encoder, the coefficient 1 < rj < 3 is specified by the the given processor [fT71 and parameter 
rpmax U pp er bounds the energy T s k that can be used by the source encoder. Function ff (D k ) 
is a classical rate-distortion function lfT6l . For the model described in this example, assuming 
that the source is independent identically distributed (i.i.d) in time with U k ,i ~ M (0, d max ) , the 




f q (Dk,T s , k 



) = jXf?(D k )xf 2 (T S!k ), 



(3) 
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rate-distortion trade-off is given by 

fi ( D k) = (log d ™ x _ ^ mmse \ , where d mmse = \ -f— + qj , (4) 

where d mmse is the estimation minimum MSE (MMSE) for the estimate of U kji given R ki lTT8ll . 
Notice that the distortion D k is upper bounded by d max and lower bounded by d mmse . 

Example 2. The sensor observes M samples of a first-order Gaussian Markov source 
[£4,i, U kj2 , — , C4,m] G K m with correlation function given by E[U k U k+j ] = d max Q k , where 
parameter < Q k < 1 is the correlation coefficient between the samples measured in slot k. 
Notice that the larger Q k is, the easier it is for the compressor to reduce the bit rate for a given 
distortion due to the increased correlation. Adapting results from @, if the source encoder uses 
a transform encoder [19], the optimal compressor produces a number of bits equal to 

P (D k7 T s>k ) = - x [f, (D k ) + f 2 (T s , fc )] + , (5) 

with fx (D k ) = log where parameter ( > 1 depends on the type of quantizer, and 

tt{T s , k )=\og(l-q 2 )* Ts % V, \ (6) 

J-s,k 

with given parameter v, which sets a lower bound on the energy T s>k as T Sjk > u/b. Equations 
©-© are obtained by assuming, similar to [6], that the energy required for source compression 
is proportional to the size of the transform encoder. Finally, notice that, since the compression 
rate must be positive, D k is upper bounded by (dmax (1 — q 2 ) Ts ' k 

B. Problem Definition 

At each time slot k, a resource manager must determine the distortion D k and the energies 
T Sjk and T t>k to be allocated to the source and channel encoder, respectively. The decision is taken 
according to a policy n := {7r fc } fc>1 , where n k := {D k ,T Syk (5* fc ) ,T t . k (S k )} determines 
parameters (D k , T s>k , T tjk ) as a function of the present and past states S k = {Si, . . . , S k } of the 
system, where the Si = {Ei, Xi,Qi, Hi} accounts for the state of the available energy Ei , for 
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the data buffer X t , for the the source observation state Qi and the channel state H^. We define 
the set of all policies as II. Policies can be optimized according to different criteria. In Sec. [Ill] 
we adopt stability under an average distortion criterion as criterion of interest, while Sec. [IV] 
addresses the optimization of the trade-off between distortion and delay. 

III. Stability Under a Distortion Constraint 

In this section, we adopt as performance criterion the stability of the data queue connecting 
source and channel encoders. We also impose the constraint that the policy guarantees the 
following condition on the long-term distortion: 



for a fixed maximum average distortion level D tolerated by the system. We define a policy as 
D -feasible if it guarantees the stability of the data queue connecting source and channel encoders 
under the average distortion constraint (fT). Recall that stability of the data queue holds if the 
distribution of X k is asymptotically stationary and proper, i.e., Pr(X fe = oo) — > |[20l . 

A. Distortion-Optimal Energy-Neutral Class of Policies 

For a given distortion D, our objective in this section is to identify a class of policies that is 
able to stabilize the data queue and satisfy the distortion constraint © as long as this is possible. 
We refer to this class of policies as distortion-optimal energy-neutral. Notice that this definition 
generalizes that of "throughput optimal" policies [1211 considered in related works such as flT], 
where only the stability constraint is imposed. By definition, a distortion-optimal energy-neutral 
class of policies n do C n contains at least one D-feasible policy. For instance, the set II of 
all policies is clearly distortion-optimal energy-neutral. However, this is a rather unsatisfying 
solution to the problem. In fact, it does not help in any way to identify a D-feasible policy for 
a given system setup. Instead, we want to identify a smaller class H do , which is parametrized 
in a way that makes it easy to evaluate a D-feasible policy. The propositions below identify a 




(7) 



fc=i 



January 18, 2013 



DRAFT 



9 

distortion-optimal energy-neutral class of policies for the separate source and channel encoders 
model depicted in Fig. \T\ and described in Sec. [n] 

Proposition 3. For a given distortion D, a necessary condition for the existence of a D-feasible 
policy is the existence of a set of parameters D q > 0, T q > for q E Q, T/ 1 > for h e H, 
and < a < 1 such that 

E Pr (?)/ 9 ( Dq i Tq s) < E Pr (%" ( T t h ) . E Fi ^° q ^ D > < 8) 

q h q 

E Pr(g)Tf < (1 - a) E [£ fc ] , and E Pr ( < " E [#*] • (9) 

(J h 

Remark 4. Parameters .D 9 , T% , T/ 1 and a, whose existence is necessary for the existence of 
a Z)-feasible policy according to Proposition |3] have a simple interpretation. In particular, Tj? 
, Z) 9 can be read as the average energy and distortion that the source encoder selects when 
the observation state is Q k = q, whereas T/ 1 can be seen as the average energy that channel 
encoder draws from the available energy for transmission when the channel state is H k = h. 
Moreover, condition ©-left is necessary for the stability of the data queue, condition ©-right is 
necessary to satisfy the constraint ©, and conditions © are necessary for energy neutrality. This 
interpretation will be used below to derive a class of distortion-optimal energy-neutral policies. 

Proof: The processes / (D k ,T S:k ,Q k ) and g(H k ,T ttk ) must be asymptotically stationary 
ergodic for queue © to be asymptotically stationary. Hence, the policy n must be asymptotically 
stationary. Under this assumption, the necessary condition for the distribution of X k to be 
asymptotically proper is [/ (D k , T Sjk , Q k )] < E„- [g (H k , T tik )} from standard results on G/G/l 
queues (see any reference on queuing theory, e.g.. [|20l Ch.3]). Notice that the average E^, with 
respect to the joint distribution of the state variables E k , Q k , H k , is explicitly dependent on the 
policy 7Tfc. From this condition, since / is separately convex in D k , T Sjk and g is concave in T tik , 
we have the following necessary condition J2 q Fr(q)f q (E^ [D k \ Q k = q] , E^ [T Sik \ Q k = q]) < 
< ^ h Pr(h)g h (E,r [T tjk \ H k = h]), where we have used Jensen inequality on both sides. Defin- 
ing D q = E,,- [D k \Q k = q], T* = E^ [T S:k \ Q k = q], and T t h = E w [T ttk \ H k = h], the condition 
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©-left is then proved. As for ® , we consider that, from ©, we must have ^ J2k=i + ^*,k) — 
~kJ2k=iEk + for if > 1, and the initial state of the energy buffer E . Then, for a 
stationary ergodic policy tt, we get E^ [T s>k ] +E n [T t)k ] < E [E k ], where ± J2 k =i T s,k -> E^ [T S)fc ], 

SfLi ?f,fc E tt [Tt, fc ], and ^ Xllli E k + jf ~^ E t^]- Given tne definitions and the inequality 
above, © are proved, having set a = E„- [T^] /E . To conclude, for ©-right, we observe 
that the distortion constraint © is satisfied. ■ 

We now look for a distortion-optimal energy-neutral class of policies. To this end, based on 
Proposition [3l it is enough to exhibit a class of policies such that it contains a D-feasible policy 
as long as the necessary conditions ®-@ are satisfied for some set of parameters D q , T q , 
T/ 1 and a. Proposition |3] suggests that it is possible to find D-feasible policies that select D k 
and T Sik based on the observation state Q k only, whereas the selection of T ttk depends on the 
channel state H k only. Based on this consideration, let us define the class of policies Tl do 



IT 



do 



D k = D q , T S)fc = min 



[l-a)E k -e,T* 



T t h = min 



aE k — e, T t 



for Q k = q 

(10) 

for H k = h 



where D q > 0, T q > for q E Q, T t h > for h E K, and < a < 1 are fixed design 
parameters. 

Proposition 5. A policy in U do is D-feasible if conditions (|#]) hold, along with 

^2Pi(q)T*<(l-a)E[E k ]-e, and Pr(g)lf < aE [E k ] — e. (11) 

Remark 6. The sufficient conditions in Proposition [5] for the policies in U do to be D-feasible 
coincide, for e — > 0, with the necessary conditions derived in Proposition |3] Therefore U do con- 
tains a D-feasible policy any time the necessary conditions of Proposition [3] hold. As discussed 
above, this implies that the set U do is a distortion-optimal energy -neutral class. Moreover, it 
should be noted that the class IT io , given (flOl) . is parametrized by a small number of parameters 
and the policies in Il do perform separate resource allocation optimizations for the source and 
channel encoders. In particular, the energy allocated to the source encoder T Sjk only depends on 
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the observation state Qk, and not on the channel quality H k , whereas the energy T t>k for the 
channel encoder only depends on H k , and not on Q k . The energy allocation between the two 
encoders is governed by a single parameter < a < 1. This entails that, once this parameter is 
fixed, and thus the energy budget available at the two encoders is fixed, resource allocation at 
the two encoders can be done separately without loss of optimality. 

Proof: For < a < 1 such that £ 9 Pr(g)T« < (1 - a) E [ E k] ~ e and E g Pr (?) T t ft < 
aE [E k ] — e, with e small, we obtain that Pr(E k = oo) = 1 asymptotically. This is true since 
E [T Sj k + T tj k] < E [Ek] in the system (0Q), so that the energy harvested is larger than the energy 
consumed on average and the energy queue is not stable [|20l Ch.3] (see also Q for the same 
argument). This leads to the asymptotically infinite size of the stored energy, as the buffer 
capacity is assumed infinite. Therefore we have, T Sjk (Q k — q) — > T q and T tjfc (H k = h) — > T/ 1 
from (fT0T > . Notice that this argument shows that in (fTOT ) one can substitute a for any number 
between and 1 leading to the same sufficient conditions ([8]) and (fTTT) . Due to the stationarity and 
ergodicity of the processes Q k and H k , f (D k ,T s>k ,Q k ) and g(H k ,T ttk ) are stationary ergodic, 
and ^2 Pr(q)f q (D q ,T q ) < Y,h^ T ( h )9 h { T t) is tne sufficient condition for the stability of the 
queue X k Ch.3]. Being lim^^ \ Y2=i D k = Pr(q)D q , for D q such that J2 q Pr(q)D q < 
D the class of policies Ii do satisfies the constraint ©. ■ 

Remark 7. A problem of interest is to find the minimal distortion D for which the set of distortion- 
optimal energy-neutral policies U do is not empty. In other words, assessing the minimal distortion 
that can be supported without causing the data queue to be unstable. Given the separate nature of 
the source and channel energy allocations, it can be seen that one should optimize both terms in 
©-left separately, once the optimal value for a has been found. In particular, when g (H k ,T tjk ) 
is the Shannon capacity, the policy T t)k that minimizes D is the water-filling [|T6l . 

B. Suboptimal Classes of Policies 

In Sec. IIII-Al a distortion-optimal energy-neutral class Ii do has been identified. This class of 
policies, as made clear by the proof of Proposition |5] requires infinite energy storage capabilities 
at the sensor node. Let us instead consider the class of greedy policies H subl that do not use the 

January 18, 2013 DRAFT 



12 



energy buffer but allocates all the energy arrival E k to source and channel coding according to 
a fraction < a q ' h < 1 that depends on both source Q k = Q and channel Hk = h states: 

USUbl { A = D q ' h , T Sjfc = ofl» h E k , T t>k = (1 - afl> h )E k for Q k = q and H k = h (12) 

where the distortion D q ' h > also depends on both source and channel states. Notice that this 
is unlike the class of distortion-optimal energy-neutral policies (fTOl) in which, as explained in 
Remark [61 energy allocation is done independently for source (only based on Q k ) and channel 
decoder (only based on H k ). Here, parameters a q ' h , D q ' h are selected on the basis of both channel 
and source states Q k and H k to partially compensate for the loss due to the greedy approach. 
For further reference, we also consider the subclass of policies W ah2 := U subl \ aq h=a f or a vj , h y 
for which the power allocation is not adapted to the channel and observation states. 

Proposition 8. Policies in Y[ subl are D -feasible if the following conditions hold: 

J2 Pr (?) Pr WE [f q {D q '\ a q > h E k )] <J2Y, Fr ^ Pr ^) E W ft 1 " ofl ' h ) E ^)] > ( 13 ) 

q h q h 

and Yl Pr fa) Px ( h ) Dq,H < D, (14) 

q h 



where the expectation E in (|T3l) is over the energy harvesting process E k . 

Remark 9. In general, the set of policies H subl is not guaranteed to be a distortion-optimal 
energy-neutral class, since the necessary conditions of Proposition [3] could hold where the 
sufficient conditions of Proposition [8] do not. This is also confirmed via numerical simula- 
tions in Sec. Illl-Dl However, for constant observation and channel states, i.e., H k = h and 
Qk = Qo f° r a U k, and for / and g linear in T s>h and T tjk , respectively, the class IP" 61 
is distortion-optimal energy-neutral. In fact, under these assumptions, the sufficient condition 
O becomes f qo (D qo ' h °, a qo ' h °E [E k ]) < g h ° ((l - a qo > h °) E [E k ]). Defining T q = a qo ' h °E [E k ], 
T t h = (1 - a qo ' ho ) E[E k ] and D q = D qo ' h °, conditions ([T3]>-([14]) correspond to ©. Thus, the 
class of policies H subl is distortion-optimal energy-neutral. 
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Proof: Due to the stationarity and ergodicity of the processes Qk and H k , the parameters 
Dk, T S) k and T tj k, are also stationary ergodic and (Q~3} is a sufficient condition for the stability of 
the queue X k M Ch.3]. Since lim^ \ YZ=i D k = E^E^?) P^)^ f ° r D<1 ' h such 
that J2 q Eh P r (?) P*(h)D q ' h < D the class of policies H subl satisfies the constraint ©. ■ 
The greedy policies introduced above do not make use of the energy buffer at all, whereas 
the distortion-optimal energy-neutral class of policies U do does. For comparison purposes, it is 
interesting to consider hybrid policies that differ from those in U do as the energy buffer is used 
only either for compression or for transmission. The first policies Yi hybl require an energy buffer 
for the channel encoder only in order to adapt the transmission power to the channel state, 
i.e., T 4j fc = mm[aEk, for Hk = h. The energy allocated to the source encoder is instead 
independent of the observation state, i.e., T s fc = (1 — a) E k . Viceversa, the second policies 
Y[h y b2 are adapted to the observation state instead of the channel state, and require an energy 
buffer only for the source encoder, i.e., T s ^ = min[(l — a) E k , T£] for Qk = q and T tt k = otEk- 

Proposition 10. Policies in n^ 61 are D-feasible if conditions ftEty-right and HTil-right hold, along 
with J2 q Pr(?)E [f (D«, (1 - a) E k )\ < J2h Pr (% h { T l 1 )- Similarly, policies n^ fe2 are D- 
feasible if conditions ftEty-right and ^U\-left hold, along with ^ g Pr(g)/ 9 (D 9 , Tf) < 
<J2 h Pr(h)E [g h (aE k )]. 

Proof: Proof follows from the proofs of Propositions |5] and |8] ■ 

C. Analog Transmission 

In this section, we consider for performance comparison an alternative class of strategies in 
which the sampled source is transmitted directly via analog modulation (see, e.g., [fT5l ). In other 
words, a block of source samples is scaled and transmitted in one slot, so as to consume T t k 
transmission energy per channel use during slot k. Energy T t ^ is selected as T tj & = mm[Ek — 
e, Tt' h ] for Qk = q and H k = h for given parameters T^' h > 0, so that it depends on the current 
source and channel states. If the bandwidth ratio b = N/M is larger than one, i.e., there are 
more channel uses than source samples, the extra N — M source samples are unused. Instead, 
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if b < 1, then a fraction 1 — b of source samples is not transmitted. Notice that this class of 
strategies does not fall in the category depicted in Fig. [Hand discussed above, since it does not 
have separate encoders. 

We assume the observation model of Example 1 with an i.i.d. source U ki i ~ M (0, d max ) and 
an AWGN channel with SNR Hk- For a bandwidth ratio of b, it is not difficult to obtain that 
the MMSE at the receiver is given by 

+ ^r-\ l for 6 > 1 



j frp s~) tt \ ) \ "^t.fcQifc+Qfc+l d max i /i c\ 

0'mmse{ 1 t,k,^k, -tik) — \ _ ± ■ U->J 

Notice that, for fairness, the average energy used for the transmission of one sample is bT t ^ 
if b > 1. Also, notice that if 6 < 1, the maximum distortion d ma:r is accrued on the fraction 
(1 — b) of samples that are not transmitted. For simplicity, we assume that analog transmission 
has negligible power spent for source acquisition, i.e., T Sj k = 0, though this is not entirely 
correct given that even in this case there is a need for sensing, sampling and analog-to-digital 
conversions. Nonetheless, these power consumption terms are also neglected in Examples 1 and 
2. Under this assumption, we have the following. 

Proposition 11. Analog transmission satisfies the distortion constraint © if the following 
conditions are satisfied: 

£ Pr (<?) PT ( h ) d mmse fa* , Q = q,H = hj<D, and £ £ Pr(g) Yi{h)T^ h < E [E k ] . 

q h q h 

(16) 

Proof: Follows similar to Proposition [6j ■ 
We can also consider a greedy policy T t k = E k , for which the power allocation is not adapted 
to the channel and observation states and energy storage is not required. 

D. Numerical Results 

In this section we compare numerically the performance of the optimal and suboptimal source- 
channel coding policies presented so far, along with analog transmission strategies. 
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Consider first a scenario where the observation and channel states are constant, i.e., Q k = q 
and H k = h for all k. The energy arrival E k has mean 1 Joule/channel use and uniform pdf 
between and 2 Joule/channel use. We consider model ® with T™ ax = 1 Joule/source sample, 
efficiency parameters ( — 1 and 7] = 1.5 for the source encoder, and the complex AWGN 
channel Shannon capacity g h (T t)k ) = N x log(l + hT t:k ). In Fig. [2l we identify the values of 
source and channel SNRs (q, h) for which different policies are able to stabilize the data queue 
and guarantee average distortion D = 0.8. We refer to these regions as "achievable regions". 
Achievable regions are given in Fig. |2]by the area above the corresponding lines. We use standard 
tools of convex optimization for their numerical evaluation. 

In Fig.[2]-b, we can further observe that the achievable regions of the distortion-optimal energy- 
neutral class (fTOl are significantly larger than those of the greedy policies (fl"2l) due to possibility 
to store energy and thus allocate resources more effectively. Moreover, by considering also the 
hybrid policies H h y bl and H hyb2 , we can see that most of the gains are obtained, in this example, 
by exploiting the energy buffer in order to allocate energy over time to the source encoder, 
whereas the gains accrued by using the battery for data transmission are less significant. This is 
observed by noticing that the achievable region of the class U do is close to that obtained by hybrid 
policies n^ 2 , but much larger than that obtained by hybrid class of policies U. hybl . The relative 
comparison between the two hybrid policies, and thus between the use of the battery for source 
or channel encoding, depends on the functions f (D k ,T Stk ,Q k ) and g (H k ,T tjk ). For instance, 
setting a lower T T s nax would change the presented results by penalizing more the strategies that 
are not using the energy buffer for channel transmission. 

We now consider the performance of analog transmission. As it is well known from rate- 
distortion theory [[T8l . for bandwidth ratio b = 1 and the given (Gaussian) source and channel 
models, analog transmission is rate-distortion optimal. Separate source-channel coding is also 
optimal (for any b > 0) if compression is assumed not to consume any energy. Here, instead, 
the achievable region of analog transmission is expected to be larger than that of strategies that 
employ separate source-channel coding, as source encoding energy costs are taken into account 
in the given model ©. This is confirmed by Fig. |2}b. However, for sufficiently larger or smaller 
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bandwidth ratios, the extra energy spent for compression is not enough to overcome the rate- 
distortion gains attained by separate source-channel coding versus analog transmission. This is 
apparent from Fig. [2]-a and Fig. |2]-c where the analog transmission is outperformed. 

We now consider a scenario where source and channel states are not constant but vary with 
two possible states, namely Q = {10 _1 , 10 2 } (Fig. [3]-a) or Q = {10~ 0,2 , 1} (Fig. [3]-b) for source 
SNR, and H = (1(T\ 10 2 ) (Fig. da) or H = (3.5, 7) (Fig. gj-b) for channel SNR. The worst- 
case observation SNR (e.g., Qk = 1CT 1 for Fig. |3]-a) and the worst-case channel SNR (e.g., 
H k = 1CT 1 for Fig. |3]-a) have probabilities p q w and p\, respectively. In Fig. |3l achievability 
regions are the sets of probability values (p^-iP^) for which different policies guarantee queue 
stability and average distortion D = 0.8. In particular, the regions are identified as the area 
below the corresponding curves. The results emphasize the importance of jointly adapting the 
resource allocation to both source and channel states in case of a greedy policy that does not 
employ the battery. This is seen by comparing the performance of the greedy schemes H subl (fT2l) . 
which adapts the policy to the current states, and n sufc2 , which does not. Moreover, comparing 
Fig. [3]-b with Fig. [3]-a, it is seen that the better "worst-case" state allows the distortion-optimal 
energy-neutral policy to satisfy stability and average distortion constraints for larger values of 
the probabilities On the contrary, the greedy policies suffer from the worse "best-case" 

state (q, h) = (1, 7) of Fig. |3]-b, as this corresponds to operating critically close to the border of 
their achievable regions (see the achievable region of H subl in Fig. |2rb). 



The stability criterion considered in Sec. |ITT] does not provide any guarantee on the delay 
experienced by the reconstruction of the source in a certain time-slot. In some applications, 
instead, one may be willing to trade distortion for a shorter delay. In this section, we address 
such requirement by looking for policies that minimize a weighted sum of distortion and delay. 
In particular, we propose to minimize the expected total discounted cost [|22i 



IV. Delay-Distortion Optimization 




(17) 



k=0 
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where < A < 1 is the discount factor and < 7 < 1. The latter parameter weights the 
importance of distortion versus delay in the optimization criterion. Notice that if 7 = then 
one minimizes the average length of the data queue, which, by Little's theorem, is the same as 
minimizing the average delay. 

In order to tackle the minimization of (fTTt over the policies n defined in Sec. [Till we 
assume that: (i) The data and energy buffers are finite; (ii) The set of possible decisions 
n k '■= {D k ,T S)k ,T ttk } is discrete; (Hi) The energy arrival E k takes values in a discrete and finite 
set; (z'v) The sets of values assumed by rates / (D k ,T S:k ,Q k ), g (H k ,T tik ), and by the queue 
length X k are discrete. Following standard theory [E2l Ch. 6], these assumptions entail that the 
optimal policy is deterministic and stationary (Markovian). In other words, {D k , T Stk , T t>k ] are 
function of the present state S k = {E k , X k ,Q k , H k } only. Therefore, the solution can be found 
via value iteration Il22l . Notice that, due to (/), data buffer overflow may happen, in which case 
the compression bits are lost and a maximum distortion d max is accrued for the current slot. 

While in general the optimal policy allocates resources to source and channel encoder through 
parameters Tf' h as a function of both source Q k = q and channel H k = h states, the class of 
policies Ii do (flOl) performs such allocation independently for source and channel encoders. For 
comparison purposes, we evaluate also the performance in terms of criterion (Q/7) of a class of 
policies that optimize separately the source encoder parameters {D k ,T Sik } as a function of Q k , 
and the channel encoder parameter {T t k } as a function of H k . As for the distortion-optimal 
energy-neutral class of policies, the energy resources are split between the encoders, such that 
the source encoder makes use of a fraction a of the energy, whereas the rest is utilized by 
the channel encoder. Specifically, the energy-buffer is divided into two buffers, that are used 
independently by the encoders: the source encoder buffer is charged by aE k , while the channel 
encoder buffer absorbs the remaining quantity of energy arrival (1 — a) E k . The source encoder 
policy {D k ,T s k } is optimized via value iteration with respect to the criterion (fTTT ) by assuming 
a constant transmission rate g (H k ,T ttk ) = g. On the other hand, the channel encoder policy 
{T t fe } is optimized via value iteration with respect to criterion ([T7T ) with weight 7 = (since 
it cannot optimize its policy with respect to the distortion), by assuming a constant source rate 
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f (D k ,T Sjk ,Q k ) — /• The best "separable" policy is finally obtained by selecting the values 
(a,g, f) that achieve the best delay-distortion trade-off. 

A. Numerical Results 

In this section we compute numerically the trade-off between delay and distortion by min- 
imizing (fTTT) for different values of 7. Specifically, for each 7, we evaluate the average delay, 
which is measured by the average data queue length by Little's law, and average distortion. 
The optimal policies are computed via value iteration [22] and so are the suboptimal policies 
corresponding to separate optimization of source and channel encoders. 

In Fig. |4] the delay-distortion trade-off is shown both for the optimal policies and for the 
"separable" ones discussed above. The discount factor is A = 0.5. The compression model 
is ©, with minimum required energy per sample v = 0.1 Joule/sample and bandwidth ratio 
6=1. The quantities of interest are discretized as follows: X k E {0, . . . , 5} is expressed in 
multiples of the codeword length M = N; The energy buffer size is E k — E k E {0, 1,2} and 
E k E {1,2} with p e w being the probability that E k — 1 (worst case); The source correlation 
values are Q = {0.1, 0.5} and channel SNR values are H = (0.5, 10), with probabilities p^ and 
for Q k = 0.1 and H k = 0.5 (worst cases); The distortion takes values as D k E {0.1, 0.55, 1} 
and d max = 1; The source-encoder rate / (D k , T Sjk , Q k ) /M is rounded to the smallest following 
integer, while the channel-encoder rate g (H k , T ttk ) /N is rounded to the largest previous integer. 

In Fig. H] we observe that the optimal policies obtain a remarkably better delay-distortion 
trade-off compared to the separable policies, both for low and large worst-case probabilities, 
i.e., Pw = Pw = Pw = Pw = 0.1 and 0.9. This demonstrates the importance of a joint resource 
allocation over the encoders whenever the delay is also of interest. Note that for increasing 
average buffer length (delay), since the buffer size is finite, it becomes more crucial to adopt 
a joint resource allocation. This is because the separate approach is not able to prevent buffer 
overflow as the source encoder operates without channel state information. 
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V. Multiple Access 



In this section, we briefly discuss an extension of the analysis to a scenario in which L 
sensors access a single access point employing TDMA. Random access protocols will be con- 
sidered in future work (an analysis with exogenous bit arrivals, and thus no source encoder, 
can be found in |fT3l ). Each sensor is modeled as described in Sec. HU and we assume that 
observation qualities Q fc = [Qi t k, • • • , Ql^Y £ Q L = [Qi, ■ ■ ■ , Ql], channel qualities H fc = 
[Hi t k, • • • , H Ltk \ T G H L = [Hi, . . . , Hi], and the energy arrivals E fc G are jointly stationary 
and ergodic. We tackle the problem of designing policies defined, extending Sec. lII-Bl as the tuple 
v := {rfc, 7^k}k>i mat consists of the scheduling policy Tfe G {1, . . . , L}, which reserves the slot 
k to one sensor I G {1, . . . , L}, and of the joint resource allocation policy 7r fc = [7r ljfc , . . . , 7r i fc ] T , 
where each entry tt/^ is defined as in Sec. III-BI and corresponds to the distortion and energy 
allocation for the I th sensor. Recall that {r^, 7r^} generally depends on the whole history of past 
and current states (see Sec. III-BI) . Note that time-slot k is exclusively assigned to sensor I, i.e., 
Tt,i,k — if Tfc 7^ I. We define a policy v as D-feasible if it guarantees the stability of all data 
queues and average distortion constraints lim^oo - Ylk=i Dk,i < Du collected for notational 
convenience in vector D = \Di, . . . , Dl\ . As for the single sensor scenario, we are interested 
in finding a distortion-optimal energy-neutral class of policies T do C T, i.e., a subset of all 
possible scheduling policies T that contains at least one D-feasible policy. 

In the following we state a necessary condition for the existence of a D-feasible policy v. 

Proposition 12. For a set of distortion constraints D, a necessary condition for the existence 
of a ID-feasible policy v is the existence of the set of parameters Df > 0, > for qi G Qi, 
T t h j > for h { G Hi, < /3/ 1 < 1 , for h G U L , and < a t < 1, such that 




(18) 



L 





(19) 



i=i 



qi&Ql 
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Pr(a)25 < (1 - a,) E [^ >fc ] , anJ ]T £ Pr(h)/3 / h I* < a,E [J5, iJfc ] . 

(20) 

Remark 13. The interpretation of Proposition [T2l is similar to the one of Proposition [3] given in 
Remark @] The additional parameters /3 z h can be interpreted as the fraction of the subset of time 
slots with joint channel state equal to h for which the sensor I is scheduled. 

Proof: A necessary condition for the stability of the queue of the I th sensor is E„ [/ (Di )k , T s j tk , Qi, 
E v [g (Hi t f., T t j } k)] Il20l . Since we have T t ^ k = for r k ^ I, on the right of the inequality we ob- 
tain E v [g (H ljk , T tjhk )] = *LheH t J2heH^Mi)=h Pr « ( T k = Z, H fc = h) x 
xE„ [g (Hi jk ,T tj i >k ) |t& = l,Hi )k = h\. Then, using Jensen inequality at both sides, the stability 
condition becomes 

P*(Qi)f qi (E. [D hk I Q hk = qi ],E v [T 8fltk | Q l)k = qi }) < E v [f {D ltk , T a , l>k , Q l>k )} < 

<J2 P l ( Tk = Z ' Hfc = h ) E - \9 ( H i,k, T t ,i,k) \r k = I, H ltk = h{\ <(21) 

<J2 Pr ( h ) F l ( Tfc = z l Hfc = h ) 9 hl ( E « [TtiAn = I, H ljk = h\). 

hi&Hi heH L Ml)=hi 

Finally, condition d) is obtained for Df = E v [D ljk \ Q hk = = E v [T Stlik | Q l>k = g,], 

/?h = Pr ^ ( Tfc = Z|H fc = h), and 7§ = E v [T ttLk \r k = l,H hk = hi]. Note that (fT9|)-left follows 
immediately from this definition. As for conditions (|20l . recall that, from (Q3, we must have 
i Ef=i (^.j,* + T tM ) < i Ef=i ^i,fc + ir, for ^ > 1, and the initial state of the energy 
buffer Ei j0 . Thus, for a stationary ergodic policy we get 

(22) 

where ± £f =1 ^V,,* ^ [T s ,^], i £f =1 £? I)fc +^ -> E [£^],and ± £)f =1 T Mjfe E„ [T tj/jfc ] = 

= E^eHi Ehe-Hi:h(0=^ Pr ( h ) Pr « ( Tfc = ^l H * = h ) E « l T t,i,k\ T k = I, Hi, k = hi}. Given the above 
definition of ffi, and T t h j, inequality becomes E 3i eQ ! Pr (*) T 5+ 
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+ E hl enX^Mi)=hM h )^ T tJ < E[J5 lffc ], proving & where a, = E„ [T u , h ] /E [J5? IfJ J. 
To conclude, condition (fT9l)-right follows from the distortion constraints similar to Proposition 

m ■ 

In order to define a distortion-optimal energy-neutral class of policies, Proposition [T2l suggests 
to consider a class of scheduling policies T do in which scheduling is done opportunistically based 
on the channel states h of all sensors according to a probability distribution ffi: if the channels 
are equal to h, then sensor I is selected with probability (3^. Notice that scheduling is independent 
of the observation qualities in a given slot. Moreover, energy and distortion allocations at each 
sensors are similar to the policies ir k in the class U do (fTOl , and, thus, in particular perform 
separate resource allocation over source and channel encoders. It can be shown that, similar to 
Proposition |5J the so defined class of policies T do is distortion-optimal energy-neutral in T. 

A. Numerical Results 

In this section we assess numerically the performance of the distortion-optimal energy-neutral 
class of scheduling policies T do . For comparison purposes, we introduce the suboptimal class 
of policies T sub , that schedules each sensor according to a fixed probability 0i = Pr (r fc = I), 
independently of the current channel conditions. We consider L = 2 sensors, which are modeled 
as in Fig. [3]-b (see Sec. IIII-DI) . For sensor 2, the probabilities (p^ 2 ,p^ 2 ) of the worst observation 
and channel states are fixed, whereas, for sensor 1, (p^,p!^) are varied. Fig. |5] shows the 
corresponding achievability region, defined, as in Sec. IIII-DI as the set of (p^,p!^) for which 
the given policy is able to stabilize the data queues and guarantee the given average distortions. 
For further comparison, Fig. |5]also shows the outer bound to the achievability region given by 
the case where only sensor 1 is present. The numerical results confirm that the achievability 
region of the optimal class of strategies T do (dashed lines) is larger than that of the suboptimal 
class T sub (dot-dashed lines). Moreover, note that the achievability regions shrink if the worst 
case states probabilities (p^ 2 ,p^ 2 ) of the second sensor get larger, since sensor 2 requires more 
transmission resources to compensate for both the worst observation and channel conditions. It 
is further interesting to observe that, for p^ = 0.1 and p^ 1 = 1, the achievability regions of T do 
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and T sub , in terms of p^, are practically the same (circle marker). This is due to the fact that 
the variations of channels are not large enough to enables gains by adapting parameters ffi. 

VI. Conclusions 

We studied energy management for a system consisting of a single sensor whose task is that 
of reporting the measure of a phenomenon to a receiver. The main problem is that of allocating 
energy between the source and the channel encoders based on the current amount of available 
energy, state of the data queue, quality of the measurement and of the wireless channel. We 
first look for a distortion-optimal energy-neutral subset of all policies, that contains at least one 
policy able to stabilize the data queue and to satisfy a maximum average distortion constraint. 
We found that optimal policies according to this criterion operate a separate energy allocation 
of source and channel encoder. Instead, we showed that a joint energy management over source 
and channel encoder is required to achieve the desired trade-off between delay and distortion. 
Finally, we considered a system with multiple sensors and obtained TDMA scheduling policies 
that guarantee the stability of all data queues, whenever the distortion constraints are feasible. 
Overall, our results, which also include further comparisons with a number of suboptimal policies, 
shed light on the challenges and design issues that characterize modern cyber-physical systems. 
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Figure 1. An energy-harvesting sensor composed of a cascade of a source and a channel encoder powered by a resource 
manager that allocates the energy available in the buffer (e.g., battery or capacitor). 
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Figure 2. Achievable regions (regions above the curves) for the digital policies (lines with markers) and for the analog policies 
(only lines): (a) bandwidth ratio b = 0.201; (b) bandwidth ratio 6 = 1; (c) bandwidth ratio 6 = 5 (Ek ~ U (0, 2), with mean 1 
Joule/channel use; D — 0.8; compression model with T™ ax = 1 Joule/source sample, £ = 1, r\ = 1.5, maximum distortion 

dmax — 1). 
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Figure 3. Achievable regions (regions below the curves) of the digital policies II do Q0>, IT"" Gil, and n"" 62 , with two 
channel and observation SNR states, respectively. (Ek ~ U (0, 2), with mean 1 Joule/channel use; D = 0.8; compression model 
{3}, with T^ ax — 1 Joule/source sample, f = 1, rj = 1.5, bandwidth ratio 6=1, maximum distortion d ma x=l)- 
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Figure 4. Delay-distortion trade-off, where average delay is proportional to the depicted average data queue length (maximum 
data-buffer length: 5 codeword lengths; maximum distortion: dmox=l> discount factor: A = 0.5; compression model ((5}, with 
minimum required energy per sample v = 0.1 Joule/sample; bandwidth ratio: 6=1; queue length: Xk £ {0, . . . , 5} expressed 
in multiples of M; energy buffer size: Ek — Ek £ {0, 1, 2}; energy arrival: Ek £ {1, 2} Joule/sample; source correlation values: 
Q = {0.1, 0.5}; channel SNR values: H = (0.5, 10); distortion values: D k £ {0.1, 0.55, 1}). 
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Figure 5. Achievable regions (regions below the curves) of the scheduling policies T do (dashed lines) and T sub (dot-dashed 
lines). The solid line corresponds to the optimal policies when the second sensor I — 2 is not present (Ei^ ~ U (0,2), with 
mean 1 Joule/channel use; Di = 0.8; compression model l[3}, with, for both sensors, T™ ax — 1 Joule/source sample, f = 1, 
r\ = 1.5, bandwidth ratio 6 = 1, maximum distortion d max =^)- 
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